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Based on rigorous scattering theory we establish a systematic methodology for research of 
metamaterial-modified current-carrying conductors, from which we mathematically demonstrate the 
explanation of transformation optics could be extended in metamaterial-modified wireless power 
transfer system, and based on that we could establish a equivalent model. More important, our 
demonstration reveals that the equivalent model will still be applicable even when TO could not 
give a direct explanation, as the requirements of complementary media is not satisfied. And numer¬ 
ical results from our methodology as well as COMSOL verified our findings. The demonstration is 
not under specific frequency, the conclusion could be extended to a broad range of wavelength, and 
expected to be applicable for active cloak etc. 

PACS numbers: 42.25.Fx, 41.20.Jb, 88.80.ht 

I. INTRODUCTION 

Metamaterials HQ are widely applied in many ar¬ 
eas @Q; in particular, they have established the experi¬ 
mental foundation for transformation optics (TO) [7H9|. 

Moreover metamaterial-modified wireless power transfer 
(wpt) Hi[13 are attractive for researchers, so as the 
development in active cloak devices EMU- 

It is well known that TO works well for the manipula¬ 
tion of EM waves ini®. Most of the time TO will 
give a direct geometric explanation, which may greatly 
simplify the explanation for metamaterial’s modification. 

While in metamaterial modified WPT, the TO is not well 
applied, as WPT system are based on current-carrying 
conducting wire and power transfer, which means those 
results based on perfect electric conductor (PEC) could 
not be applied directly in WPT systems. Rather than 
numerical analysis from finite element method (FEM), 
to establish rigorous TO explanation in such system, an 
analytical analysis is essential. What’s more, in WPT or 
active cloak researches, the objects (e.g. receiver) might 
destroy the condition the TO explanation will rely on. 

As a result it is hard to extend TO explanation in such 
situation. 

As a matter of fact, we develop a general method¬ 
ology based on rigorous scattering theory to ad¬ 
dress metamaterial-modified current-carrying conduc¬ 
tors, from which we are able to discuss the TO explana¬ 
tion mathematically. We choose super-scatterer [p| as 
the metamedia to modify the WPT. To clarify our con¬ 
clusions we discuss a well simplified WPT system, with 
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only two parallel conducting wires acting as an receiver 
and emitter, we call it PCW in this article. Then we cover 
the emitter wire with super scatterer shell (ss-shell), we 
name it ss-PCW. With introduced applied voltage, we are 
able to discuss its power transfer properties. Based on 
this heuristic model and our methodology, we could es¬ 
tablish an analytical equivalent model, which is extended 
from TO. It will greatly simplify the analysis of the ss- 
PCW, as it could avoid solving for detail EM fields in 
metamaterial, when we calculate the transfer power and 
efficiency of this system. Moreover, from our demonstra¬ 
tion we find the equivalent model could still be applicable 
even when the requirements of complementary media [l8] 
is not satisfied. Which greatly expand the usage scope of 
the equivalent model. 

We started from discussing metamaterial modified 
current-carrying conductor because it is a heuristic model 
for WPT and active cloak. Although the results are ob¬ 
tained from the two-dimensional (2D) circular cylindri¬ 
cal case, it is expected to be applicable for a more gen¬ 
eral cases, even three-dimensional (3D) models. And our 
analytical demonstration is not under any certain fre¬ 
quency, the conclusion could be applied to a broad range 
of frequencies. Thus the deduction from our methodology 
could be easily extended to many other research fields. 

This article is organized as follows: First, we con¬ 
sider the EM radiation and scattering problems for a 
super-scatterer shell (ss-shell)-covered conducting wire, 
in which current is introduced by impedance boundary 
conditions, where TO could be well extended. Second, 
we study a more complicated model, the ss-PCW. From 
which we give analytical analysis to discuss the limitation 
and extension of TO, and further analysis reveal that the 
equivalent model is still applicable even when TO could 
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not be applied directly. Finally, we verify our demonstra¬ 
tion through numerical results obtained from FEM and 
our series expansion solution. 


II. TO EXPLANATION FOR SS-SHELL 
COVERED CONDUCTING WIRE 

To present TO explanation on an metamaterial modi¬ 
fied wireless power transfer system, we begin with a sim¬ 
ple but heuristic model of an ss-shell-covered conductor, 
as shown in FigureUJa). It is a 2D model, and we divide 
it into three sections: section 1 (0 < r < n), section 2 
(ri < r < 7 * 2 ) and section 3 ( 7*2 < r < 7 * 3 ).In cylindrical 
coordinates, as shown in Refs. @,0 , each section has 
different permittivity and permeability tensors, 

( £ - u ~ 

l £0 = ne = jhif'ir), ( 1 ) 

{ £z= Hz = 

which can be derived from the following coordinate trans¬ 
formations: 

{ /i(r) e 7/r, 0 < r < ri, 

f 2 (r) = nn + (r ~rmf ( rm~ r2) , r 1 <r<r 2 , 

hir) = r, r 2 <r < r 3 . 

(2) 

where m/0. 

From the viewpoint of TO (complementary media), 
section 3 (vacuum) will be complementary to section 2, 
which is an ss-shell filled with negative index material 
(NIM), and we will ultimately detect an amplified image 
originating from section 1 in the domain 0 < r < 7 * 3 . 

This result can also be proven using the scattering the¬ 
ory of EM waves. As demonstrated in Ref. 0 , if we 
consider only a transverse electric (TE)-polarized EM 
field with a harmonic time dependence of e~ luJt in our 
2D model, the general series solution for the electric field 
(and the magnetic field) can be written as 

00 

E z (r,e)= J2 [anJn(k 0 f(r))+b n H^(k 0 f(r))]e ind , 

n= — 00 

(3) 

( 2 ) 

where fco is the wave vector in vacuum and J n and Hn J 
are the nth-order Bessel function and the nth-order Han- 
kel function of the second kind, respectively. 

As mentioned above, we can consider that any EM 
boundary inside the ss-shell is amplified from the do¬ 
main 0 < r < 7*1 to 0 < r < 7 * 3 , with an amplification 
factor of rj = If we place a conducting wire of radius 
R in section 1, it will form an enlarged image with an 
amplified boundary in the domain 0 < r < 7 * 3 , as shown 
in FigureHJWe typically consider a PEC surface when 
designing illusionary optical devices; here, however, we 
consider the general case of a current-carrying conduct¬ 
ing wire (to which a time-varying voltage is applied). Be¬ 
fore we solve for the parameters of the series expansion 



FIG. 1. Simple schematic of a metamaterial-modified current- 
carrying conducting wire: (a),ss-shell covered conducting 
wire. (b),The conducting wire and ss-shell are not concen¬ 
tric. 


solution, we should apply impedance boundaries to this 
model. As discussed in Ref. [19| - the conductor we con¬ 
sider here has high conductivity, and the wavelength is 
also sufficiently large to satisfy the impedance boundary 
approximation. In fact, it is also not essential to solve 
the Maxwell equations inside the conductor to obtain the 
exact solution for the EM fields in that region; as a result, 
applying impedance boundary conditions could simplify 
our calculations. 

Using the Maxwell equations and impedance bound¬ 
ary conditions, we can write the boundary condition for 
boundary 1 (the interface between the conducting wire 
and section 1 in Figure a) ) as follows, 


m n (r) 


E c (r) - E c 

^zn V J sz 


\r=R 



(4a) 


H 0n( r ) 


-iuj/j.0 1 , 


(4b) 


where n = 0, ±1, ±2,..., and the parameters £ c , // c , cr c , 
and E c szn denote the material permittivity, permeability, 
and conductivity and the induced electric held (induced 
by the applied AC power source), respectively. Moreover, 
fie 1 denotes the 6 components of the permeability tensors 
of the materials in section 1, which can be evaluated from 
equations CD- 

In each section, the series expansion solution for the 
EM fields can be expressed as 


00 

£“(r»= y C(r)e^ 

n= — 00 


00 

= E K^n(fco/a(r))+ 6 “i/i 2 ) (fc 0 /a(r))]e^,( 5 a) 

n = —00 


tf?(r,*)= E HZ n (r)e in t 


with Hg n (r)= ^ 2) Wa(r))] <(5b) 


IUJ flQ a 
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where the upper or lower index a is equal to 1,2 or and Hq associated with section 1, 2 or 3, respectively. In 
3 and denotes the corresponding / a (r), /i0 a , a", b", this manner, we can obtain the following three boundary 

equations for boundary 1 to boundary 3: 

I 


KUkoMr)) + b\H^\koh{r))} \ r = R = E c zn (R), 

/ 1 '^)klJ:(fco/i(’-))+fejg< 2) '(fco/i(r))] 

—iujfiei 


r = R = -(E c zn (R) - E %J 


gc-J- 


(6a) 


K4(^o/i(r)) + b l n H^\k 0 h{r))} | r=ri 

/i , (r)[^.J , :(fco/i(r))+^g„ ( 2 ) '(fco/i(r))] | 

-iuntn l r=r i 


[«n Jn(k 0 f 2 (r)) + b 2 n H^ ] (k 0 f 2 {r))\ 
^(’•)[^J:(fcofe(’-))+i'^n (2), (io/2(r))] 


—icone 2 


|r=ri 5 
|r=n 5 


(6b) 


[a^ J n (fc 0 / 2 ( r )) + b 2 n H ( ?\k 0 f 2 {r))} 

/2( r )[ Q n J n( fe o/2(^))+^ ij ~n (2)/ (feo/2(r))] 


r=r2 — 
r=r2 = 


where n = 0,±l,±2,.... In theory, we can obtain all of 
the EM fields by solving the above equations obtained 
from the boundary conditions. 

It is meaningful to study the EM fields for r = r]R , 
considering that this position is simply the amplified 
boundary of the conducting wire from TO theory, as 
shown in FigureOJa). Here, we deduce from equation 
© that /i(.R) = h(riR) = r]R, fi(n) = /2(ri) = r 3 and 
^2(^2) = = ^2, all of which are easy to demon¬ 
strate. As a result, we have equations ^( r ) = , 

and aforementioned equations (l6aj) , (l6b|) and (l6c|) above 
can be simplified as follows: 


i fr 1 = b 2 = b 3 

t u n u n u n j 

where n = 0, ±1, ±2 ,.... 

Moreover, by comparing equation 
©, we ultimately obtain 


(7) 


aj) with equations 


Eln(riR)= [a 3 n J n (k 0 f 3 (r )) + b 3 n H&(k 0 f 3 (r))] \ r=vR 


E c zn {R\ 


(8a) 


hL(vR)= 


f 3 ( r )[ a lJ' n (k 0 f 3 {r)) + b 3 n H n {2) ' (k 0 f 3 (r))] 


-ILd/lQS 


~{E c zn (R) - Et 


V 


yr—qR 


(8b) 


where n = 0,±1,±2,.... These equations imply that 
there is an impedance boundary at r = r]R , 


g?» | 

( 4 M - |r_r '* 

Indeed, if we take 




(9) 


o' — T] 2 <7 + lUJ£ c (rj 2 — 1) 


( 10 ) 


[<4 Jn{koh{r)) + fo3^ 2 )(fc 0 / 3 ( r ))] | r=r2 , 

/3( r )[ Q n J ^( fe o/3(^)) + ^^n (2)/ (feo/3(r))] I 

—iune 3 k—r 2 5 


equation ([9|) will have the same form as equation (Hal) . In 
other words, we can obtain an enlarged current-carrying 
conducting wire, as indicated by TO theory. 

The equivalence is discussed in detail for the following 
two cases. If there are observers at r > r]R , they will de¬ 
tect the same EM fields from the following two emitters: 

a) A conducting wire with radius i?, conductivity a 
and applied voltage E c szn that is covered by an ss- 
shell as shown in Figure Ufa). 

b) A conducting wire with radius r]R , the conductiv¬ 
ity g' defined in equation m and applied voltage 

JZ/ szn % 

Based on this equivalence, we refer to case b as the ma¬ 
terialized model of case a later in this article. 

All of the deductions above can also be used to explain 
the model depicted in FigurelTJb). However, it should 
be remarked that the TO method may be not applicable 
under certain circumstances, as we will discuss later in 
this article. 


III. LIMITATION AND EXTENSION OF TO IN 
SS-PCW 


A. EM field in ss-PCW 

Now that we have introduced TO explanation in the ss- 
shell-covered conducting wire model from rigorous scat¬ 
tering theory, we will discuss the ss-PCW, schematically 
illustrated in Figure|2fa). 

The model can also be represented as a 3D model. An 
ss-shell-covered emitter wire forms part of a circuit (so 
are the receiver wire) with an time-varying applied volt¬ 
age, which is equivalent to introduce a time-varying elec¬ 
tromotive force by applying an external EM filed to the 
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FIG. 2. Schematic illustrations for ss-PCW:(a), a simplified 
WPT system with its emitter covered by an ss-shell, receiver 
could be located in any positions out of section 2 and 1. 
(b),materialized model for ss-PCW. (c) translation of wave 
functions in different coordinates. 


emitter circuit [20j. When electric energy is transferred 
from the emitter to the receiver, it can be regarded as a 
radiation and scattering problem of EM waves. For con¬ 
venience in calculation, all conducting wires are circular 
cylinders; thus, that our previously applied methodology 
can also be used here as a suitable tool for calculating 
the transfer power and efficiency. 

In the calculation, different receiver positions require 
different series expansion solutions. Unlike Figure JTJ we 
divide section 3 into two parts: section 3a and section 3b. 
Both of them are in vacuum; section 3a is the domain 
with its radius smaller than the minimal radius which 
could cover the amplified emitter boundary. We will dis¬ 
cuss the difference for receiver wire inside section 3a or 
outside of section 3a later. 


Let us place the emitter and receiver wires in sections 
1 and 3, respectively. As in equation 0 above, we can 
write the EM fields in sections 1 to 3 as follows: 

oo 

E?\pA)= E WuMk 0 h{\~t\))e in ^ ( 11 a) 

n= — oo 

+blH^\k 0 h{\-f -p5Zt\))e in + BX ] 

OO 

4 2) (p,0)= ]T [d 2 n Mk 0 f 2 (\tm (lib) 

n = —oo 

+g*H<?\koM\?\))]e in +, 

OO 

E?\p,<t>) = Y, {alH£\k 0 h(\~t\))e in + ( 11 c) 

n= — oo 

+b 3 n Hi 2 \kof 3 (\t ~ P^\))e in4,B2 ], 

Here, we present only the expressions for the electric 
field; the corresponding Hq components can be derived 
using the equation Hq = — . 

It is important to emphasize that the electric field (in 
sections 1 and 3) consists of three parts, each of which 
has a series expansion solution in different cylindrical co¬ 
ordinates. As shown in figures [2fa) and[2fb), point O is 
the center of the ss-shell, whereas in sections 1 and 3, Ba 
is the center points of the emitter or receiver wires (where 
a = 1,2). In section 1 or 3, the origin points O, and Ba 
offer two different cylindrical coordinate systems. 

Following Weng Cho Chew’s book [2l| , we present the 
translation of the wave functions corresponding to origin 
points A and B as follows (shown in Figure j2fc)): 

Hm\ko\pl\)e irn ^ b = H£\k 0 \pt - p^\)e im ^ 

E H^l m {k 0 ^ b )J n {k 0 ^a)e in ^-< n - m )^\ ( P a<Pab) 

n =~™ ( 12 ) 

E Jn-m(k 0 ^ b )H£\k 0 rt)e in ^-*( n - m ^, ( Pa >Pab) 

n = —oo 

where m = 0, ±1, ±2,.... These equations describe the 
translation from origin point B to origin point A, and 
similar equations can be used for translation between any 
two of the three points, and similar for J n . To solve for 
the expansion parameters of these series solutions, we 
need to apply the EM-field boundary conditions at each 
boundary. 

Based on the translation equations given above, as pre¬ 
viously discussed, we can write the boundary conditions 
as (with n, m = 0, ±1, ±2,...) 


MkovRi) E [aUm-nikovro^e-^-^] + I&H™ (kovRi) = E^{Ri), 


m= — oo 


JtikovRi) E [am^-n(fet)roi)e“ i(n “ mWl ] + biH^'(k oV RQ = f - E„ n (R{)), 


(13a) 


m= —oo 


Jn{k 0 R 2 ) E KH^_ n (kor 0 2)e-^ n - m ^ 2 } + b 3 n H^ 2 \k 0 R 2 ) = E^(R 2 ), 

m= — oo 

JniW £ KH^koro^e-^-^} + b 3 n H^' (k 0 R 2 ) = ico^^E^(R 2 ), 


(13b) 
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{ d^Jmikorjn) + g^ n H^\korjn) 
^m^m(^oryri) + (. k 0 rjn ) 


a^J m (fc 0 77 n) + H^\korjn) [bnJm-n(korjroi)e z(m n) ^ >1 ], 

n= — oo 

am^m(^oryri) + (fco^i) X) [^n^m-nC^oryroOe - ^ 771-71 ^ 1 ], 

n= — oo 


di l J m (k 0 r 2 ) + fhrMrn (k<ir2 ) 
d 2 m J'm(kar 2 ) (fcor 2 ) 


a^^ 2) (fc 0 r 2 ) +JmCfcora) E [&ME(to»'o 2 )e~ <(m -” ) * , ]> 

n= —oo 

a 3 m H£ y (k 0 r 2 ) + J' m (k 0 r 2 ) E [&^E Jfcoro^e-^-™^ 2 ]. 


(13c) 


(13d) 


Here, we treat the emitter as conducting wire and set 
its center at points B1 in section 1; thus, we have 

^oi = ||pom|^K0om)|Toi = IpobiI? 0i = 0om- Con¬ 
ducting wire has radius i?i, conductivity <ti and an n-th 
order applied voltage of E szn { in this article, this param¬ 
eter appears as the voltage per meter in the 2D model); 
moreover, we add a variable E^n to denote the electric 
field close to an emitter wire. For the receiver in section 

3 (or outside section 3), we have the corresponding con- 

( 2 ) 

stants do 2 , ^ 02 ,02, R 2 , and the a 2 and the variable E zn • 
Besides, the time-varying voltage applied to the conduct¬ 
ing wire of the emitter is E szn ; we set a non-zero value 
for this voltage only at zeroth order, thereby obtaining 
E szn = E sz So^ n . Similar as deduction of equation © we 
will have 


[ flm= E [b 3 n H^_ n (koro2)e-^ m - n ^}, 

J n= —00 

\ 00 

\ a 3 m = E KJm-n(koVroi)e-^ m - n ^], 

\ n= — 00 


(14) 


where ra = 0, ±1, ±2,-For each order n, we can sim¬ 

plify these eight equations into four equations and thus 
have only four unknown coefficients. It should be em¬ 
phasized the deduction m could be applied for receiver 
in section 3a and 3b or out of section 3. 


B. Materialized model for ss-PCW 

We have previously mentioned that the materialized 
model for ss-shell covered conducting wire. By applying 
TO method we could also provide the materialized model 
here for ss-PCW. The schematic of ss-PCW’s material¬ 
ized model is shown in Figure. ©b), the electric field in 
the materialized model consists of two components: 

00 

E z {PaAa)= E b 3 m H ^(k 0 \-f-^2 |)e im ^ 2 

m= — oo 
00 

+ E blH^ikotf -mt\)e in ^(15) 

n= — oo 

These two components are expressed in two different co¬ 
ordinate systems with different origin points correspond¬ 
ing to the center points of each wire. 

However, it should be emphasized that when the re¬ 
ceiver wire is placed in section 3, the TO explanation 
will lose its utility because the domain (section 3) used 
for the complementary ss-shell (section 2) is not in vac¬ 
uum (i.e., the conditions for complementary media are 
not satisfied). This phenomenon has also been noted 
in Ref. [l5]. Fortunately, under certain specific circum¬ 
stances, our materialized model will still be applicable 
even when the complementary media requirements do not 
hold. We demonstrate this as follows. 

When the receiver wire is outside of section 3a (ro 2 = 
Pob 2 > Pob x = roi), we can write the emitter wire’s 
boundary conditions as 


Jn(koR'l) E [aLJrn-nMje-^-™^] + 


m =—oo 


JUkoR'i) E [«m Jm-n(k 0 r' 01 )e- i(n - m) ^} + b 1 n Hi 2 ) '(k 0 R[) = iuJ - ^d^ (El n (R[) - E szn (R[)), 


m= — oo 


pi 


where a} m = E l b k H m-k( k o r 02 )e l(m fe) E (r 02 > r' 01 ), 


(16a) 

(16b) 

(16c) 


k= — oo 


where n = 0, ±1, ±2,..., and the parameters are R[ 


r]Ri (the radius of the emitter wire in the materialized 
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FIG. 3. comparison for ss-PCW and its materialized model: 
(a), when receiver are out of section 3, complementary me¬ 
dia is not destroyed, the TO will give a direct explanation for 
the equivalence between these two models (the conductivity of 
materialized emitter in materialized model should be trans¬ 
lated as we deduced), (b), receiver situated in section 3b, 
the complementary media is destroyed, however our demon¬ 
stration indicate that the equivalence still exist, (c), receiver 
penetrated in section 3b, the equivalence is invalid. 


model), p 0 B 1 = r' 01 = i?r 0 i < r 02 , and <!>[ = <j>\. Sim¬ 
ilar to the deductions presented above, we set o' as in 
equation m, and choose the same E szn as in equations 
(|13ap . Finally, we find that equations (I16a|) and (I16b[) 
are identical to (toil) and m- We should note that the 
relations of equations (j!6c|) are satisfied given the condi¬ 
tion ro 2 > 7*01. When ro 2 < r f Q1 , as shown in equations 
CD, the translation for equations (flbl) will be different. 
Which also means that the receiver should be outside of 
section 3b when we consider the equivalent materialized 
model. 


C. Limitation and Extension of TO 

In conclusion, when the receiver wire penetrates into 
section 3b (as shown in Figurej2ja)), the TO explanation 
(complementary media) will not function correctly, but 
the equations for the coefficients in ss-PCW model will 
remain equivalent to its materialized model (shown in 
Figure |2fb)), which also means that equations (H6j ) are 
still applicable. The relationship between TO, ss-PCW 
and it’s materialized model are presented in Figurej3] 

It is interesting to consider the case in which the re¬ 
ceiver wire is moved even closer to the emitter, for exam¬ 
ple, penetrating into section 3a or the amplified image of 
the emitter wire’s boundary. Where, equation will 

not be satisfied, and thus, there is no simple geometric 


explanation as our materialized model based on TO. To 
calculate the EM fields, we must solve equations (fl3l) of 
emitter and those of receiver numerically. 


D. verify the conclusion numerically 

Analytical demonstration is presented above, here we 
will verify our conclusion numerically through solving 
equation CC3> or applying FEM in COMSOL. Once we 
obtain the EM fields in the model, the total power in¬ 
troduced to the system, transfer power and transfer effi¬ 
ciency are calculated as follows 

^total = MIKA - E sz )2nRj £ -^Et-±-}, 

(17a) 

^receiver = E Re ^E 2 z * n E 2 zn 2nR 2 

n= — oo 

(17b) 

c = P] g ceiver . (17c) 

Motal 

As presented in the above demonstration, we have al¬ 
ready demonstrate the TO based materialized could be 
equivalent to its ss-PCW even when section 2 is not 
complementary to section 3. As the equivalence also 
lead to the equivalent transfer power and efficiency, we 
could verify our demonstration numerically by using the 
FEM (COMSOL) and our series expansion solution. We 
present the numerical results in two different scenarios, 
in which receiver is situated in section 3b and section 3a 
respectively. 

First, we consider the case in which the receiver wire 
is located in section 3b, as shown in Figurel^b), which 
means that 77 ( 7*01 + Ri) < ro 2 - Here we calculate the 
model for roi = 0.1 m, cfoi = 0.1 m, Ri = 0.01 m, 
R 2 = 0.03 m, ^02 = 0.5 m, and ro 2 = 0.96 m, and 
set the amplification factor to 77 = 6 (we call this sce¬ 
nario ss-PCW 1). Because 77 ( 7*01 + Ri) < 7 * 02 , we can 
consider a equivalent materialized model with only re¬ 
ceiver and emitter wire with parameters of r' 01 = 0.6 m, 
d' 01 = 0.6 m, R[ = 0.06 m, R 2 = 0.03 m, d 02 = 0.5 m 
and ro 2 = 0.96 m. The voltage applied to the emitter 
wire is 10 V/m (and the wavelength is 100 m). The ma¬ 
terial of the emitter wire is conducting (with a = 5 x 10 5 
S/m), and the receiver wire is connected to an external 
load (0.6 + 0.6zfl/m). The results are listed in TableHl As 
we expected, although the TO approach loses its validity 
when the receiver wire enter section 3b, the material¬ 
ized model remains useful. Figurej4ja),(b),(c) shows the 
electric fields (z component) in the three simulated mod¬ 
els. By analyzing the distribution of E z , we will find the 
equivalence between ss-PCW and its materialized model 
easily. 
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TABLE I. Results from COMSOL and the series expansion solution for ss-PCW 1 




Series expansion 

solution 

COMSOL results 


Power (W/m) 

Ratio(%) 

Power (W/m) 

Ratio(%) 


Receiver 

Emitter 

Efficiency 

Receiver 

Emitter 

Efficiency 

ss-PCW 

0.227 

0.559 

40.6% 

0.223 

0.538 

41.4% 

materialized model 

0.227 

0.559 

40.6% 

0.225 

0.543 

41.4% 

ordinary PCW 

0.057 

0.229 

24.7% 

0.056 

0.221 

25.3% 


TABLE II. Results from COMSOL and the series expansion solution for ss-PCW 2 




Series expansion 

solution 

COMSOL results 


Power (W/m) 

Ratio(%) 

Power {W/m) 

Ratio(%) 


Receiver 

Emitter 

Efficiency 

Receiver 

Emitter 

Efficiency 

ss-PCW 

3.82 

4.39 

87.0% 

3.67 

4.27 

86 . 0 % 

materialized model 

0.368 

0.853 

43.1% 

0.368 

0.851 

43 . 2 % 

ordinary PCW 

0.021 

0.235 

9.1% 

0.021 

0.235 

9 . 1 % 



FIG. 4. Electric fields (z component) for ss-PCW 1 & 2: 
(a),electric fields (z component) in ss-PCW 1. (b),electric 

fields (z component) in materialized model of ss-PCW 1. 

(c) ,ordinary PCW. Compare (a) with (b), we will find the 
equivalence between ss-PCW 1 and its materialized model. 

(d) ,electric fields (z component) in ss-PCW 2. (e),electric 

fields (z component) in materialized model of ss-PCW 2. 
(f),electric fields (z component) in ordinary PCW. Compare 
(d) with (e), we will find there isn’t equivalence between ss- 
PCW 2 and its materialized model. 


Similar as the above presented comparison, when re¬ 
ceiver wire penetrated into section 3a (we call it ss-PCW 
2 , in fact ss-PCW 2 can not be represented by its re¬ 
lated materialized model as we demonstrated analyti¬ 
cally above). Based on our series expansion solution and 
COMSOL results we could also give the comparison as 
Table HI In the ss-PCW 2 scenario, we set the parame¬ 
ters as roi = 0.5 m, do i = 0.3 m, R\ = 0.05, R 2 = 0.2 
m, do 2 = 0.5 m and ro 2 = 2.49 m, with amplifi¬ 
cation factor as 77 = 6 , which satisfies the condition 
77 ( 7^01 + -R1) > 7“ 02 • The wire’s material and external load 
are taken to be the same as in ss-PCW 1 , so as the ap- 


|| 9645 plied voltage.Figure J4||d),(e),(f) shows the electric fields 
4 (z component) for ss-PCW 2 scenario. 

2 In conclusion, it is obvious that the scenario ss-PCW 

2 cannot be explained in terms of the TO approach and 
also cannot be correctly represented by a corresponding 
materialized model, while ss-PCW 1 shows the TO based 
materialized model still be applicable even the condition 
of complementary media is destroyed. 

Moreover, it should be noted, because of non¬ 
monotonic transforms, the computational error in sim¬ 
ulation for superscatterers from FEM is hard to control, 
this phenomenon is also discussed in Ref.[22[. To give 
the above mentioned meaningful results the computa¬ 
tional time for FEM is much longer than that cost by 
our methodology. Besides, FEM requires an increasing 
number of mesh elements as the scale of the model grows, 
the time required for compuation and the iterative error 
will increase even if there is sufficient memory. Our re¬ 
sults obtained from series expansion solution are more 
reliable for simulation of super scatterer, and need less 
computational time, especially for large-scale model. 

IV. CONCLUSION AND DISCUSSION 

In this article we extend the TO explanation in meta¬ 
material modified WPT system, and from rigorous scat¬ 
tering analysis, we establish a equivalent model for the 
system, named materialized model. From our analytical 
analysis, we demonstrate that the TO based materialized 
model could be applied even when TO approach could 
not give a direct explanation, where the requirements for 
complementary are not satisfied. Besides, we present the 
numerical results from FEM and our series expansion so¬ 
lution to verify our findings. The conclusion is firmly 
verified by these numerical results. 

The wireless power transfer system here is simplified 
as two parallel conductor wires, however, the results and 
methodology should also be useful for a more complicated 
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model, such as more pairs of wires or coils, and even a 
3D helix coil. 

Moreover, our methodology which is based on multiple 
scattering theory could deal with the calculation of EM 
fields in superscatterer modified WPT more efficiently 
than FEM, especially in large-scale model. 

Furthermore, our analytical demonstration is not un¬ 
der some certain frequency, which means the results could 
be applied in researches under a broad range of wave¬ 
length, e.g. the research in active cloak and antenna. 
The demonstration was under impedance boundary con¬ 
dition, which could be applied in many researches. We 


expect the analysing for a more specific boundary condi¬ 
tion like perfect electric conductor boundary could reveal 
more special properties by applying our methodology. 
Besides, the TO explanation for WPT is very important, 
which might significantly improve the WPT research. 
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